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We study the spectra and damping of surface plasmon-polaritons in double graphene layer struc- 
tures. It is shown that application of bias voltage between layers shifts the edge of plasmon ab- 
sorption associated with the interband transitions. This effect could be used in efficient plasmonic 
modulators. We reveal the influence of spatial dispersion of conductivity on plasmonic spectra and 
show that it results in the shift of cutoff frequency to the higher values. 



I. INTRODUCTION 

Recently the possibility of fabrication of double-layer 
graphene structures was reported^. In such structures 
two graphene layers separated by thin dielectric are clad 
between the two sections of substrate. A number of de- 
vices based on double-layer graphene has already been 
fabricated, among them there are field-effect tunneling 
transistors 2 and optical modulators 3 . Plasmonic prop- 
erties of double-layer graphene are of particular interest 
too: first, such structures can act as plasmonic modula- 
tors^; second, plasma resonances can affect the charac- 
teristics of optical modulators, resulting in greater mod- 
ulation depth near plasmon peak^. However, the spectra 
of plasma oscillations in double graphene in the presence 
of bias voltage between layers have not been studied yet. 

In this paper we derive the spectra of surface plasmon- 
polaritons (SPP) in double-graphene layer (2GL) struc- 
tures using the model of graphene conductivity presented 
ir£. In our calculations the spatial dispersion of conduc- 
tivity is taken into account, which turns out of great im- 
portance at the elevated Fermi energy fi comparable of 
higher than the temperature T. We analyze the effect of 
bias voltage applied between layers on the spectra and 
damping of SPP and show that the 2GL structures can 
be used as efficient plasmonic modulators. 



II. ELECTRODYNAMICS OF SURFACE 
PLASMON-POLARITONS IN 2GL 

The 2GL structure under consideration is schemati- 
cally shown in Fig. Q] It consists of two GLs separated by 
a barrier layer. One of the edges of each GL is connected 
with a contact, while the other edge is isolated. The bias 
voltage Vb is applied between the contacts. Thus, one GL 
serves as the gate for another. The dielectric constants 
of the top, middle, and bottom media are denoted by 
St, £ m , and £;,, respectively; the distance between GLs 
is d. The x-axis coincides with the direction of the SPP 
propagation parallel to the metal electrodes, the 2-axis 
is orthogonal to the graphene plane. 

Due to the peculiarities of graphene conductivity, the 



FIG. 1. Schematic view of the 2GL structure and directions 
of vectors E, H and k in the TM SPP 



propagation of TE SPP modes along its surface is pro- 
hibited^. For this reason we focus on TM-modes with 
the ac electric and magnetic fields E = {E x , 0, E z }, 
H = {0, Hy, 0}. We search for the solutions of Maxwell 
equations in the form: 



E, 



E t = E ot e l(fcr - tJt) -' ttZ , (1) 
0mt e K ™ z + E 0mi e- K ™ z ) e «(««-u.*) ) (2 ) 
E h = Eo b e i(kx ~ ut)+KbZ , (3) 



where the indices t, m, and b, again, stand for the 
top, middle and bottom media, respectively, k% = 
yjk 2 — eiLu 2 /c 2 , and c is the speed of light. 

Matching the vectors E and H on the graphene sur- 
faces (the conductivities of top and bottom layers are at 
and (76 ), one arrives at the dispersion law 
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Eq. Q is simplified for symmetric structures with e t = 
£b — £ and <T t = (Tb = a. The latter condition is valid 
provided the residual doping of graphene layers is weak. 
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In the situation at any finite bias voltage the modules of 
Fermi energies in GLs are equal but have opposite signs, 
hence, due to electron-hole symmetry the conductivities 
of layers are equal. The simplified dispersion laws are 

(^_4™V— W— ( 5 ) 

\ K J K rn \ 2 J 

( - ina ] tanh (*2*) + ^ = 0. (6) 

Eq. ([S]) corresponds to a "symmetric mode" with the cur- 
rent oscillations in both layers in phase, while Eq. ([6]) rep- 



resents an "antisymmetric mode". We are particularly 
interested in symmetric oscillations^. At small distances 
between layers {n m d <C 1) their dispersion depends only 
on the net conductivity of two layers^: 



K t K b u 

III. MODEL OF GRAPHENE CONDUCTIVITY 

To perform a numerical analysis of SPP spectra we use 
the model of graphene conductivity, presented in Rcf. 6 : 
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(8) 



Here the indices 1 and 2 stand for conductance and va- 
lence bands, in particular ei(p) = pvf and £2(p) = 
—PVf, vf — 10 6 m/s, p± = p±7ik/2, /(e) is the electron 
distribution function (the equilibrium Fermi function is 
assumed), v x = vfCOS0 p and v\2 — ivpsiii9 p are the 
matrix elements of velocity operator. The first line of 
Eq. ((5J) corresponds to the intraband transitions, while 
the second one accounts for the intraband. To account 
for the momentum relaxation the frequency should be 
treated as ui —> ui + ir^ 1 . The expressions for t~ x for 
different scattering mechanisms could be found in 10 . Be- 
ing interested in high-quality samples with /i » T wc 
restrict ourselves to the consideration of electron-phonon 
scattering. The appropriate relaxation time isi 2 - 

-1 _ -lP"F 

where t^ 1 ~ 3 • 10 11 is the characteristic electron-phonon 
collision frequency at room temperature T = 26 mcV. 



We stress that in the following calculations the spa- 
tial dispersion of conductivity could be neither omitted 9 
nor treated in a perturbative manner—. At elevated 
Fermi energies and at high dielectric constants of sur- 
rounding materials the dispersion k(uj) of SPP is close to 
the threshold line k = ui/vp, hence, the electronic tran- 
sitions are substantially indirect. The latter leads to a 
pronounced reduction in interband absorption compared 
to the direct case (the imaginary part of conductivity 
chanes too). To clarify the effect formally we pass to 
the elliptic coordinates in Eq. (J5J): 2p x — k cosh u cos v, 
2p y = fcsinhusinw. Note that the wave vector is com- 
plex and the functions |p ± Kk/2\ should be treated as 
absolute values of complex vectors. Extracting the real 
part of conductivity with the Sokhotski theorem in the 
limit r _1 — > we arrive at 
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(10) 



where the dimensionless measure of spatial dispersion q = IV. RESULTS AND DISCUSSION 

\k\vF/oJ is introduced. 

The spectra Refc(w) and Imk(u>) of SPP calculated 
with Eqs. ([S]) and ([5]) are presented in Figs. [5] and |3] 
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FIG. 2. Dispersions Re fe(w) of SPP in 2GL at different Fermi 
energies. Inset: low-energy "photonic" behavior of spectra 
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FIG. 3. Damping Imfc(u) of SPP in 2GL at different Fermi 
energies. Inset: "quality factor" of SPP defined as Re kj Im k 
at the same Fermi energies 



The structure parameters are St = e m — = 5 (boron 
nitride 14 ), d ~ 2 nm. The three characteristic regions 
can be singled out in the spectra. 

Small wave vectors and frequencies k <C 
8(e 2 /h~c)p,/(h~c): the dispersion is purely "photonic", i.e. 
the dispersion law reads 



(11) 



Intermediate wave vectors and frequencies 
8(e 2 /hc)n/(hc) < k < w/vf, w < 2/i/7i. In this 
situation the coupling with photons is negligible and 
the spatial dispersion of conductivity is negligible too. 
For this reason we can set K ~ k and obtain the SPP 
dispersion in the form 
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(12) 



which yields a well-known square dependence 1 ^ Re k ex 
lo 2 until the interband transitions start playing significant 
role. 
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FIG. 4. Comparison of SPP dispersions obtained with and 
without accounting for the spatial dispersion of conductivity 



High frequencies u ~ 2p/h. A simple analysis based 
on Pauli exclusion principle states that the direct inter- 
band transitions are allowed at frequencies u > 2(1/1%. 
At such frequencies the incident light is effectively ab- 
sorbed, which governs the operation of graphene-based 
optical modulators^. The SPP at such high frequencies 
are also expected to undergo strong damping. However, 
the situation for SPP is different as the interband transi- 
tions in the case are substantially indirect. From Eq. (|9]) 
it is readily seen that the interband absorption at finite 
k is not as strong as at k = 0; in the threshold case 
|k| > u/vf the interband transitions are forbidden at all 
due to nonconservation of energy and momentum. On 
the contrary, intraband transitions at |k| > u/vp are ex- 
tremely strong as Keaintra exhibits a singular behavior 

Rea in tra OC [(kv F /u)) 2 - Vf 1 ! 2 . 

The lack of interband transitions at finite wave vec- 
tors leads to the higher cutoff frequencies of SPP in 
graphene compared to photons. The situation is illus- 
trated in Fig. HI where the two approaches for the cal- 
culation of plasmonic spectra are compared. The dashed 
line is plotted without accounting for the spatial disper- 
sion of conductivity, while in the solid curve it is taken 
into account. It should be empathized that even the in- 
trinsic spatial dispersion due to nonzero wave vector of 
SPP significantly affects the spectra. An extrinsic spa- 
tial dispersion of conductivity (e.g. created by spatially 
periodic gates 1 ^ or surface roughness 1 ^) can also change 
the characteristics of SPP cardinally. A simple analysis 
of relation (IT21 yields that the higher the permittivity 
e and the Fermi energy p are, the more important the 
intrinsic spatial dispersion is. 

Despite the shift of SPP cutoff frequency to higher val- 
ues, the modulation of propagation length by applying 
bias voltage Vb between layers is still efficient. To study 
this dependence we implement the plate-capacitor model 
relating Vb to the charge density p in layers 



And 



{V b - 2/i) = p; 



(13) 
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FIG. 5. Illustration of SPP propagation with frequency 
uj/2n = 25 THz in 2GL structure (d = 2 nm) at different 
bias voltages. The colors correspond to the amplitude E x . 

positive voltage VJ, corresponds to the positive charge 
density. In the limit /j>T the latter is given by 



(14) 



The Eqs. (flllHl~4l allow us to evaluate the damping of 
SPP of given frequency as a function of bias voltage. The 
dependence is illustrated in Fig. [5] both the propagation 
length of SPP and "quality factor" Refc/Imfc decrease 
with decreasing the bias voltage. At almost zero bias the 
wave damping could be even greater than it is predicted 



by our equations as the electron-hole scattering will play 
significant roleii. 

The plasmon modulators based on 2GL are techno- 
logically convenient as they do not require metal gates 4 . 
Moreover, the plasmonic spectra in gated graphene struc- 
tures are linear with velocity much lower than speed of 
lights, thus coupling of plasmons and electromagnetic 
radiation could not be efficient. The double layer struc- 
tures have no such limitations. 



V. CONCLUSIONS 

We have obtained the spectra of surface plasmon- 
polaritons in double-layer graphene structures. The 
waves can exhibit high " quality factors" Re kj Im k or- 
der of 10 2 at Fermi energies above 150 meV. The damp- 
ing of the waves is mostly due to interband transitions, 
which can be tuned by applying bias voltage between lay- 
ers. Those interband transitions are strongly indirect, i.e. 
accounting for spatial dispersion of conductivity is cru- 
cial for describing the behavior of SPP. We have shown 
that spatial dispersion of conductivity shifts the cutoff 
frequency of SPP to higher values compared to 2[i/h ex- 
pected from analysis of Pauli blocking. 
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